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PREFACE 

This  technical  report  has  been  prepared  specifically  in 
connection  with  the  Seminar  on  Theoretical  Wave-Resistance  to  be 
held  at  the  University  of  Michigan  during  August  1963.  This  re¬ 
port  contains  both  new  theoretical  results,  previously  unreported, 
and  material  previously  presented  in  HYDRONAUTICS  Technical  Re¬ 
port  117-2. 
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NOTATION 


A  (0),A  (0) 
1  2 

A(0),B(0) 


B 

f,f 


g 

H 

k 

o 

k 

i 

L 

m(x) 
m(x, z ) 

R 

S  (x),S  (x) 
1  2 

V 


X 

1 


z 

1 


Coefficients  of  polynomial  representing  source 
distribution 


Wave  amplitude  functions 

Coefficients  of  polynomial  representing  doublet 
distribution 

Beam 

Depths  of  end  points  of  doublet  distribution 

Froude  numbers  with  respect  to  draft  and  length 
respectively 

Acceleration  of  gravity 

Draft  of  ship 

=  gH/V2 

=  gL/V2 

Length  of  ship 

Total  strength  of  source  per  unit  length 
Total  strength  of  source  per  unit  area 
Ship  wave  resistance 

Functions  of  x  defined  by  Equations  [4]  and  [5] 
Uniform  velocity  at  <« 

Right  handed  rectangular  coordinate  system  with 
z  positive  upward,  x  in  the  direction  of  the 
uniform  flow  velocity  V,  and  origin  on  the  mean 
free  surface 

Dummy  variable  indicating  x  coordinate  of  source 
distribution 

Dummy  variable  indicating  negative  z  coordinate 
of  doublet  distribution 
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€ 

n 

s ’  ^sB’ ^  s  s 
^BF 


- 


A 


Perturbation  coefficients  of  b 

n 

Wave  heights  due  to  total  ship,  bow  only,  and 
stern  only,  respectively 

Wave  heights  due  to  infinite  bulb  and  finite 
bulb  respectively 
_  r  _  f 

SB  ^BF 

Strength  of  doublet  per  unit  length 


Strength  of  quadrupole 
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INTRODUCTION 

There  are  two  ways  of  looking  at  the  merits  of  a  bulb  on  a 
ship.  One  is  its  effect  on  the  interactions  between  the  waves 
themselves  and  the  other  is  its  effect  on  the  ship  resistance  it¬ 
self  as  represented  by  these  waves.  Recently  Inui,  Takahei,  and 
Kumano  (i960)  studied  the  bulbous  bow  ship  from  the  former  point 
of  view  by  considering  the  bulb  as  a  point  doublet.  They  paid 
special  attention  to  the  l8o°  phase-difference  between  the  bulb 
wave  and  the  ship  bow  wave,  and  tried  to  make  both  amplitude  func¬ 
tions  the  same  in  order  to  cancel  out  the  wave.  Wigley  (1936) 
mentioned  this  wave  phenomenon,  but  attacked  the  problem  via  the 
wave  resistance  approach.  He  derived  six  rules  for  the  design  of 
a  bulb  for  a  ship  bow 

This  latter  approach,  involving  the  determination  of  the  wave 
resistance  directly,  has  been  the  conventional  one  used  by  most 
investigators.  However,  this  generally  leads  to  many  cumbersome 
mathematical  calculations.  On  the  other  hand,  the  method  dealing 
with  the  interaction  between  the  waves  themselves  can  be  simple 
and  leads  to  a  better  physical  insight.  Using  this  approach, 
Takahei  (i960)  studied  the  "waveless  bow"  further  and  Kumano  (i960) 
investigated  the  "waveless  stern."  They  tried  to  match  up  the 
amplitude  function  of  a  doublet,  which  is  a  function  of  the  doub¬ 
let  strength  and  the  Froude  number,  as  closely  as  possible  with 
the  amplitude  function  of  a  ship.  However,  it  was  not  very  clear 
how  close  the  doublet  distribution  could  be  matched.  It  will  be 
shown  in  this  report  that  a  doublet  distribution  whose  amplitude 
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function  is  exactly  the  same  as  that  of  the  ship  can  be  found  but 
that  it  extends  along  a  line  from  the  free  surface  to  infinite 
depth.  With  such  a  doublet  distribution  at  the  bow  and  stern 
the  waves  from  the  bow  and  stern  of  the  given  symmetric  ship  may 
be  completely  canceled  and  the  ideal  fluid  wave  resistance  of  the 
ship  becomes  zero.  The  semi  infinite  nature  of  the  extension  of 
the  doublet  distribution  need  not  be  a  serious  restriction  since 
it  will  be  shown  that  the  effect  of  discarding  that  part  of  the 
doublet  distribution  extending  from  some  finite  depth  to  infinity 
can  be  made  very  small. 

Further  it  is  shown  that  a  quadrupole  distribution  at  the 
bow  and  the  stern  in  addition  to  the  doublet  distribution  offers 
practical  advantages. 

WAVES  DUE  TO  A  SHIP  AND  A  BULB 

As  in  the  usual  analysis,  an  invlscid,  incompressible,  and 
homogeneous  fluid  in  steady  flow  wit',  a  free  surface  and  of  in¬ 
finite  depth  is  considered.  The  coordinate  system  C-xyz  is  right 
handed  with  origin  on  the  mean  free  surface,  x  positive  in  the 
direction  of  the  uniform  flow  velocity  V,  and  z  positive  in  the 
upward  direction.  In  this  report,  except  in  the  sections  on  two- 
dimensional  streamlines  the  quantities  x,  y  and  other  lengths  ex¬ 
cept  z  are  non-dimensionaii zed  with  respect  to  the  ship  length  L; 
z  is  non-dimensionalized  with  respect  to  the  draft  H,  and  m  is 
non-dlmensionalized  with  respect  to  V  and  LH.  All  equations  are 
expressed  in  non-dimensional  form.  All  wave  height  components 
are  assumed  to  be  small  and  additive. 
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It  is  well  known  that  a  point  source  of  strength  m  located 

at  a  point  (x  ,0, -z  ),  where  z  >0,  produces  a  wave  height  non- 

11  i 

dimensionalized  with  respect  to  L  at  large  x  given  by  (see  Have¬ 
lock,  1951), 

f  T/2 

£  =  8k  /  m  exp(-k  z  sec2©)  sec3©  cos  |k  (x-x  )sec© 

s  0  Qt  [ 1  i  . 

0 

X  cos  (k  y  sin  0  sec2©)  d©  [l] 

i 


where 


k  =  ^  and  k  =  ^  . 

*  V*  0  V2 


Hence,  for  a  line  distribution  of  sources  at  y  =  0, 

z  =  -z  ,  0  <  x  <1  represented  by  the  series, 

i  i 


00 


n=0 


(2] 


the  wave  height  will  be 
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tt/2  1 

C  =  Bk  /  I  m(x  )  exp(-k  z  sec2©)  sec3©  cos  (k  (x-x  )sec©) 

SOjJj.  0l  xi 

o  o 


X  cos  (k  y  sin  ©  sec2©)  dx  d© 

i  i 


8k 


i  o 


tt/2 


exp(-kQz  sec2©) 


^  sin  (k 


x  sec  0) S  (0) 


-  cos(k  x  sec  ©)S  (0)  -  sin  (k  (x-l)  sec  ©  )  S  (1) 

12  1  1 


+  cos  (  k  (x-l)  sec  ©  )  S  (l) 

1  2 


cos  (k  y  sin  ©  sec2©)  d© 


[31 


with 


OO 

-I 


(-1 )nm(^n)  (a ) 
(k  sec© )2  ^  ^ 


[4] 


n=0 


S  (a) 


■I 


,  «  \ n  (2n  +  1 ) ,  » 
(-1)  mv  '  (a ) 

/.  „»2n-l 

(k  sec©) 

1 


[5] 
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where 


m^n  ^ (a  )  = 


j  8nm(x 

v  n 
ox  l 


x  =  a 

i 


If  we  put 


^sB  +  ^ss  ’ 


where 


asB 


8k 


i  o 


tt/2 


exp(-kQZ  sec2©) 


sin(k  x  sec  0)S  (0) 

i  i 


cos(k  x  sec  0 )S  (0) 

1  2 


w/2 


cos(k  y  sin  0  sec20)  d0,  [6] 


8k  f  r 

C  =  — -  exp(-k  z  sec20)  -sin  (k  (x-l)sec  0)  S 

33  k  3  J  01  L  1  1 


(1) 


+  cos 


(k  (x-l)  sec  0 )S  (l)  cos 

1  2 


(k  y  sin  0  sec20 

i 


d0, 

[7] 


then  C,  can  be  interpreted  as  the  wave  from  the  bow  and  C  as 
S  s  s 

the  wave  from  the  stern  according  to  the  idea  of  the  elementary 

wave  (see  Havelock,  193i+a)  •  We  can  see  that  both  £  and  £  are, 

St)  s  s 
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in  general,  composed  of  sine  and  cosine  elementary  waves.  Hence¬ 
forth  we  will  omit  the  word  "elementary"  except  to  avoid  ambigu¬ 
ities.  In  the  integrand  of  £  _  and  £  ,  the  coefficient  of 

S.D  ss 

sin 

cos  (k  x  sec  0)  cos  {k  y  sin  0  sec20)  is  called  an  amplitude 
function. 

For  a  ship  of  a  finite  draft,  4  has  only  to  be  integrated 

s 

from  z  =  0  to  z  =1.  In  this  case  S, (a)  may  be  a  function  of 
ii  1 

2  . 

i 

We  now  consider  a  doublet  of  strength  -p.  nondimensionali2ed 

with  respect  to  V  and  LH2  located  at  x  =  0,  y  =  0,  z  =  -2  .  The 

i 

nondimensional  wave  height  due  to  this  doublet,  at  large  x,  is 
(see  Wigley,  1936) 

ir/2 

4B  =  -8k q2 J  M-  exp(-kQz  sec20)  sec40  sin(k  x  sec  0) 
o 

X  cos  (k  y  sin  0  sec20)  d0  [8] 

i 

where  4  is  nondimensionalized  with  respect  to  L.  Note  that  the 

sine  waves  originating  from  the  origin  have  a  negative  sign  in 
this  case. 


CONDITION  OF  NO  WAVE  RESISTANCE 

We  have  seen  from  the  previous  section  that  the  wave  height 
due  to  a  submerged  body  or  a  surface  ship  may  be  represented  at 
large  x  by 
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V2 


C  = 

o ' 


A  (0)  sin  (k  x  sec  0)  +  A  (0)  cos  (k  x  sec  0) 

XI  2  1 


X  cos  (k  y  sin  0  sec2  0)  d0 


[91 


Havelock  (193^) showed  that  the  nondimensionalized  wave  resist¬ 
ance  represented  by  this  wave  system  is 


tt/2 


R  = 


A  2(0)  +  A  2(0)  \  cos30  d0 

1  a  | 


[10] 


where  R  is  related  to  the  wave  resistance  R  by  R 

o 


The  integrand  of  this  wave  resistance  Integral  is  always 
positive.  Hence,  if  there  exists  a  ship  whose  wave  resistance 
is  zero,  the  necessary  and  sufficient  condition  for  the  no  wave 
drag  ship  is 


A  (0)  =  A  (0)  =  0,  [11] 

1  2 

(i.e.  no  wave  itself).  Of  course  the  trivial  solution  for  this 
is  the  case  of  no  singularity  or  of  an  infinitely  thin  frictionless 
flat  plate  parallel  to  the  uniform  stream,  i  e.  u  =  0  and  m  =  0. 
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However,  there  may  be  a  non-trivial  solution  which  makes 

A  (0)  =  A  (0)  =  0  with  the  proper  selection  of  the  singularity 
1  2 

distribution.  Krein  (1955)  proved  that  there  does  not  exist 
any  finite  ship  whose  Michell's  resistance  (the  wave  resistance 
represented  by  Michell's  integral)  becomes  zero,  but  he  cited  a 
few  examples  of  infinite  convoys  whose  wave  resistance  is  zero. 

In  the  next  section  we  will  find  waveless  solutions  involving 
the  appropriate  combination  of  a  doublet  distribution  with  a 
source  distribution.  The  basis  for  these  is  that  while  a  posi¬ 
tive  sine  wave  usually  starts  from  the  bow  of  a  ship  without  a 
bulb,  a  negative  sine  wave  starts  from  the  surface  point  directly 
above  a  bulb  representable  by  a  point  doublet.  This  has  been 
checked  experimentally  as  well  as  theoretically  by  many  investi¬ 
gators  . 

COMBINATION  OF  DOUBLET  AND  SOURCE  DISTRIBUTION 

The  problem  is  therefore,  how  to  find  a  doublet  distribution 
for  a  given  or  otherwise  determined  proper  source  distribution  in 
order  to  make  the  total  wave  amplitude  function  zero;  i.e.  find 

the  doublet  distribution  \i(x  ,  z  )  on  the  longitudinal  center 

i  i 

plane  such  that  for  a  given  source  distribution  m(x  , z  )  on  the 

l  i 

same  plane 


n(x  z  ) 

1  X 


—  00  —00 


exp(-k  z  sec2@)sec30  cos 

o  i 


k  (x-x 

L  1  1 


)sec  0 


dx  dz 

l  i 


exp(-k  z  sec20 )sec40sin 

o  i  7 


k  (x-x  )sec0 

i  i 


dx  dz  =  0 

i  l 


8 
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Let  us  consider  the  simplest  case  first;  vis.,  a  uniform 

source  distribution  on  the  line  0  <  x  ,  y  =  0,  z  =  f.  To  a 

11  i 

first  approximation  this  represents  a  semi-infinite  paraboloid 
and  is  equivalent  to  the  case  of  n  =  0  in  Equation  [2].  Then 
the  contribution  to  the  wave  from  the  bow  is  from  Equation  [6] 

7t/2 

8ko  f 

C  _  =  — ; —  a  exp(-k  f  sec20)sec20  sin  (k  x  sec  0) 

^sB  k  I  o  o  i  '  '  j,  ' 

o 

X  cos  (k  y  sin  0  sec20)  d0  [13] 

i 

Now  this  can  be  matched  by  a  uniform  doublet  distribution  b  on 

o 

a  vertical  line  x  =  0,  y  =  0,  f<z  <  Then  the  corresponding 

i 

wave  height  will  be,  by  an  integration  of  [8]  with  respect  to  zx 
from  f  to  00  } 


f/2 

I  b= 


=  -  8kQ  j  bQ  exp(-kQf  sec20)sec20  sin  (k  x  sec  0) 


X  cos  (k  y  sin  0  sec20)  d0 


[14] 


Prom  Equations  [13]  and  [14]  we  now  see  that 


<  =  ^SB  +  CB  =  ° 
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Thus  according  to  this  result  the  waves  from  the  point  on 
the  surface  above  the  nose  of  a  submerged  semi-infinite  para¬ 
boloid  will  be  completely  cancelled  by  a  vertical  circular  cylin¬ 
der  extending  from  the  nose  of  the  paraboloid  to  infinite  depth. 
This  is  the  fundamental  idea  of  how  the  wave  amplitude  function 
can  be  made  equal  to  zero. 

Let  us  now  consider  a  little  more  general  case  of  Equation 
[6]j  e.g.,  an  axially  symmetric  body.  However,  the  source  dis¬ 
tribution  of  Equation  [2],  which  is  an  ordinary  power  series,  is 
not  the  proper  one  to  use  for  the  doublet  distribution  because 
this  will  produce  cosine  terms  according  to  Equations  [6]  and  [7] 
as  well  as  sine  terms,  while  our  expression  for  the  waves  due  to 
the  doublet  distribution  at  x  =  C  gives  only  sine  terms.  Since 
there  is  no  easy  way  to  make  the  amplitude  function  of  the  cosine 

terms  A  (6)  of  Equations  [9]  and  [10]  equal  to  zero  by  combination 
2 

with  a  doublet  distribution,  one  thing  to  do  is  to  employ  source 
distributions  which  do  not  give  cosine  terms.  The  other  idea  is 
to  consider  a  different  singularity  distribution  which  produces 
cosine  waves  as  well  as  sine  waves  to  cancel  ship  waves  This 
will  be  discussed  later. 

With  reference  to  Equations  [2],  [8],  [5],  and  [6],  we  note 
that  m^°^(0)  =  a  m^^O)  =  2!  a - n/2  ^(0)  =  (2n)!a  and 

m^1 ^ (0)  =  a  . m^2n  +  1^(o)  =  (2n  +  l)!  a  [15] 

1  2n  + 
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i.e.,  S  (0)  Is  related  to  only  the  even  powers  of  the  series  [2] 

1 

while  S  (o)  is  related  to  only  the  odd  powers  of  the  series.  The 

2 

amplitude  function  of  the  sine  terms  in  Equation  [6]  is 


8k 

_ o 

k  3 


exp(-kQz  sec20)  S  (0), 


and  that  of  the  cosine  terms. 


exp(-k  z  sec2©)  S  (0). 
o  2 

We  now  see  that  a  source  distribution  given  by  an  even 
power  series  does  not  produce  cosine  terms.  According  to  Weier- 
strass'  approximation  theorem,  any  continuous  curve  can  be  ap¬ 
proximated  by  a  polynomial  in  the  domain  0  %  x  <  1-  If  we  con¬ 
sider  another  curve  in  -1  <  x  <  0,  symmetric  with  respect  to 
the  first  about  x  =  0  the  polynomial  which  represents  the  whole 
curve  in  -1  <  x  ^  1  must  be  an  even  power  series.  Hence  we  may 
say  any  curve  in  0  <  x  <  1  can  be  approximated  by  a  polynomial 
of  even  powers  in  the  domain.  Equation  [2]  may  therefore  be 
written 


8k 
_ o 

k  3 


n  =  0 


[16] 


11 


HYDRONAUTICS ,  Incorporated 


If  we  assume  that  our  ship  is  symmetric  about  its  midsection  at 
x  =  C.5  then  m(x  )  is  odd  with  respect  to  this  point  or 

i 

m(0)  =  -m(l),  m'(O)  =  +  m'(l) 

m(n)(0)  =  (-l)n+1  m(n>(l)  [17] 

Since  from  Equation  [16]  m^2n+1^(0)  =  0,  from  Equation  [17] 

m('"n+1)(l)  =  C.  Hence  from  Equation  [5]  S  (C)  =  S  (1)  =  C. 

2  2 

Similarly  from  Equations  [17]  and  [4]  we  can  easily  see 

S  (G)  =  -  S  (1) 

l  l 

Then  from  Equations  [6]  and  [17],  we  see  that  the  amplitude 

functions  of  i  .  and  L  are  exactly  the  same  and  both  are  sine 
Sd  s  s 

waves . 

When  the  body  is  not  symmetric  fore  and  aft  with  respect  to 

the  midsection,  the  wave  from  the  stern  can  be  treated  separately 

in  a  similar  manner  by  changing  coordinates.  If  we  take  the 

origin  at  the  stern,  and  express  the  source  distribution  for  the 

body  by  an  even  power  series,  the  wave  system  from  the  stern  can 

be  analyzed  in  exactly  the  same  way  as  for  the  bow. 

Therefore  the  solution  to  the  zero  wave  drag  problem  for 

f  can  be  treated  exactly  in  the  same  manner  as  that  for  i 
S3  SB 
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We  now  consider  only  £  of  Equation  [6]  corresponding  to 

SB 

the  source  distribution  of  Equation  [16]  along  z  =  f  ,  y  =  0, 

i  i 


8k 


tt/2 


ssB 


* 


exp(-kQf  sec2©) 


(-1)  (2n)!a 


n.o(V°e 


en..  N  sin  (k  x  sec  9) 


a (n-1 ) 


X  cos  (k  y  sin  ©  sec2© )  d©  [l8] 

i 

The  terms  for  large  n's  become  very  small  for  the  usual 
range  of  Proude  numbers,  and  in  practice  we  have  only  to  consider 
Just  a  few  terms  of  the  series  given  by  Equation  [l6]  (see  ap¬ 
pendix).  In  order  to  find  the  corresponding  doublet  distribution 
which  cancels  the  wave  height  C  ni  let  \i(z  ),  the  doublet  strength 
per  unit  length  be  expressed  as 


.(O  =  ^  bn(z^-f)n  =  bnqn  [19] 

n  =  0  n  =  0 


on  the  line  x  =  y 


0,  f  <  z  <  oo. 

i 


Thus  for  f  <  z  <  f  , 

i  l 
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ir/2  f  -f 


^BF  "  8ko2 


bn^n  sec4* 


O  0 


n=0 


exp[  -k  (q+f  )sec20 


X  sin  (kj_  x  sec  0)  cos  (k  y  sin  0  sec2©)  dq  d0 

i 


tt/2 


=  -8 


b  n! 
n 


n=0 


k 

o 


exp(-kQf  sec20)-  exp(-kQf  sec20) 


‘I 

m=f) 


sin(k  x  sec  0)cos(k  y  sin  0  sec20)d0 


[ 20, a] 


The  value  of  the  second  term  for  f  -*  °°  vanishes  (see  Appendix). 

Thus,  for  f -*  oo, 

x 


tt/2 

t_  =  -8  /  exp(-k  f  sec2©)  ,  - - - — t-t- 

■B  I  o  '/  n-l.ee*  (»-l)e 


n !  b 


sin  (k  x  sec  0) 


n=0 


X  cos  (k  y  sin  0  sec20)  d0 


[20, b) 


If  we  compare  Equations  [lb]  and  [20,b]  we  obtain 
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n 


(.if  (2n)! 
v  '  n! 


2n+l 


2n 


[21] 


as  the  condition  that  the  regular  waves  given  by  Equation  [lb] 
be  exactly  cancelled  out  by  those  given  by  Equation  [2C,b].  The 
doublet  distributions  [19]  located  at  x  =  C  and  extending  from 

z  =  f  to  z  =»  generate  a  wave  system  which  for  an  ideal  fluid 

l  i 

cancells  completely  the  bow  wave  system  due  to  any  given  source 
distribution  [Equation  16].  In  a  similar  manner  as  mentioned 
before,  the  stern  wave  system  can  be  completely  eliminated 

For  the  surface  ship,  we  have  only  to  consider,  m(x  )  of 

i 

Equation  [l6]  as  the  '.ource  strength  per  unit  area  on  the  center 
plane  y  =  C,  for  t  <  ;;  £  1,  C  <  x  <  1.  If  we  integrate  both 

l 

[lo]  and  [2C]  with  respect  to  f  from  f  =  0  to  f  =  1,  we  obtain 
the  corresponding  amplitude  functions 


i  o 


exp(-k_)f  sec20) 


I 


(-l)n(2n):a 


2n 


(k  sec  9)2 

n=u  x 


(n-l) 


[22, a] 


and 


B(e)  =  8 


exp(-k^f  sec2©) 


I- 


nib 


n=0  o 


1  2 (n- 1  )fl 

sec  9 


df 


[22, b] 
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If  Equation  [21]  holds  for  all  depths  f  we  then  get 

A(e)  -  B(e)  =  o. 


Thus  for  any  vertical  distribution  of  sources.  Equations  [19] 
and  [21]  will  give  the  required  bow  and  stern  bulb  shape  to  com¬ 
pletely  cancel  the  ship  waves. 

It  should  be  emphasized  that  the  b^  in  Equation  [22]  is  no 

longer  exactly  the  same  as  that  obtained  in  the  expansion  of  the 
doublet  p(z  )  in  Equation  [19]-  Prom  Equations  [22],  [19]  and 

i 

[20]  the  doublet  strength  per  unit  depth  is  now 


in  the  domain  0  <  z  <1,  and 

l 


>  [23, a] 


^(z  ) 

l 


C)n  dC  for  z  £  1.  J 

l 


For  the  case  in  which  the  source  strength  does  not  vary  with 
depth 


b  =  const 
n 


and 
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or 


and 


dp.(z  ) 

_ i_ 

dz 

i 


CO 


n=0 


M-(z 


“  ,  n+1 

Eb  z 

n+1 

n=0 


for  0  <  z  <1 

i 


M-(2^  ) 


n=0 


>[23, b] 


for  z  ^  1 

i 


Equations  [23,b]  show  that  the  slope  of  the  doublet  distribution 

along  the  line  x  =  0,  y  =  0  has  a  discontinuity  at  z  =1. 

x 

Especially  if  we  consider  the  one  term  n  =  0  (which  is  an  im¬ 
portant  term  for  the  practical  range  of  Froude  numbers  [see 
Appendix] )  the  bulb  shape  may  be  similar  to  a  round  pointed 
pencil . 

In  summary,  we  have  shown  how  A  (0)  and  A  (0)  of  Equation 

X  2 

[9]  can  be  made  identically  zero.  Thus  a  means  of  satisfying  the 
necessary  and  sufficient  condition  for  a  zero  wave  drag  ship  in 
an  ideal  fluid  has  been  obtained.  Unfortunately  this  requires 
the  use  of  bulbs  of  infinite  draft.  In  a  later  section  we  will 
consider  the  effect  of  cutting  off  the  doublet  distribution  (or 
bulb)  at  some  finite  depth. 
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QUADRUPOLE  DISTRIBUTIONS 

We  have  represented  the  ship  shape  in  the  preceding  by  an 
infinite  even  power  series  (l6).  However,  the  series  may  not 
converge  too  rapidly,  and  for  the  higher  Proud e  numbers,  which 
we  may  wish  to  consider,  the  higher  order  terms  in  the  series 
(l6),  can  not  be  neglected.  Suppose  we  have  an  arbitrarily 
specified  polynomial  of  finite  terms  which  represents  the  source 
distribution  on  the  center  plane  (see  Weinblum,  195t).  Then  we 
have  to  consider  the  odd  power  terms  in  Equation  [2]  and  thus 
the  cosine  wave  system  occurs  as  well  as  the  sine  wave  system  in 
Equation  [6]  and  [7]»  starting  at  the  bow  and  the  stern.  We 
therefore  have  to  consider  some  means  of  cancelling  cosine  waves. 

A  point  doublet  is  a  one  step  higher  order  singularity  than 
a  point  source;  the  wave  or  the  flow  field  due  to  a  doublet  is 
therefore  represented  by  a  derivative  of  that  due  to  a  source 
with  respect  to  the  position  of  the  source  (same  as  the  position 
of  the  doublet)  in  the  direction  of  the  doublet  (see  Milne- 
Thomson).  It  Is  well  known  that  a  submerged  point  source  pro¬ 
duces  positive  cosine  waves  (see  Equation  [l]).  Therefore,  the 
negative  sine  wave  can  be  produced  by  a  negative  doublet  which 
produces  a  closed  body  In  the  uniform  stream. 

Now  we  consider  a  one  step  higher  order  singularity  than  a 
doublet  which  is  called  a  quadrupole.  That  is,  we  consider  two 
doublets  of  opposite  signs,  with  the  magnitude  of  each  strength 
u,  and  the  distance  between  them, a,  nondimensionalized  with  re¬ 
spect  to  H.  Along  the  same  idea  of  forming  a  doublet  by  a  source 
and  a  sink  (see  Ml Ine-Thomson  1955),  we  form  a  quadrupole  by 
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making 


lim  ap.  =  finite  constant,  say,  X 
a  -*•  0 


[24] 


which  is  the  strength  of  the  quadrupole.  Its  waves  will  be  the 
derivative  of  the  wave  height  due  to  a  point  doublet  with  re¬ 
spect  to  the  position  of  the  doublet  in  the  direction  of  the 
quadrupole,  namely  cosine  waves.  The  sign  of  the  waves  will  de¬ 
pend  on  the  sign  (direction)  of  the  quadrupole.  That  is,  a 

point  quadrupole  with  the  strength  X  at  x  =  0,  y  =  0,  z  =  -z 

o  i 

in  the  uniform  flow  V  generates  the  wave  height 


X  exp(-k  z  sec20)sec50  cos(k  x  sec  0) 

o  1  1 


X  cos  (k  y  sin  0  sec20)  d0 

l 


[25] 


where  X  =  X^^L  ) 

There  are  two  great  features  of  the  quadrupole  which  will 
be  mentioned  in  the  following  sections.  One  is  the  same  feature 
which  the  doublet  distribution  had  with  respect  to  the  even 
power  terms  of  the  ship  source  distribution,  but  in  this  case 
with  respect  to  the  odd  terms.  The  other  is  the  composition  of 
a  point  singularity  which  produces  no  wave  -  say,  a  no-wave- 
singularity. 
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COMBINATION  OF  QUADRUPOLE  AND  SOURCE  DISTRIBUTIONS 

We  have  mentioned  that  the  odd  power  terms  of  the  ship 
source  distribution  [2]  produce  cosine  waves  and  that  the  quad- 
rupole  also  produces  cosine  waves.  Otherwise,  both  amplitude 
functions  have  similar  qualities  to  those  of  the  even  power 
terms  of  Expression  [2]  and  for  the  doublet.  Namely,  if  we  con¬ 
sider  the  odd  power  terms  of  Expression  [2]  say 

m  =  ax  2n+1  in  0  <  x  <  1,  [26] 

2n+l  1  1 

0  <  z  <  f, 

1 

then  the  wave  due  to  this  source  distribution  starting  from  the 
bow  x  =  0  Is  from  Equation  [6] 

tt/2 

C  =  -  8  j  a2n+1l l-exp(-kQsec2e  )] (-1  )n  (2n+l)!  X 

QJ 

X  k  2  (n+-*-  ^cos£?n+J  0  oos(k  x  sec  0)  cos(k  y  sin  0  sec20)d0 

1 

[27] 

In  exactly  similar  manner  as  in  matching  the  doublet  distribution 
to  the  even  power  terms  of  the  ship  source  distribution  in  order 
to  cancel  the  bow  waves,  the  corresponding  quadrupole  distri¬ 
bution  may  be  put  as  in  Equation  [23,b] 
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A 


n+1 


n+2 

A _ 

n+2 


in  0  S  z  <  1  [28] 

1 


r  n+2  /  Nn+2, 

[z  -  (z  -1)  ] 

X  -  b  n  = - -  in  1  5  z  <» 

n+1  n+2  1 


[29] 


then  from  Equation  [25],  the  wave  due  to  this  quadrupole  becomes 


tt/s 


n-1 


£  =  -  8k  b  n  1-exp (-k  secz0)  (n+l)i  k  X 

^q  1  /  n+1  '  o  '  '  '  o 


X  cos  y  cos(k  x  sec0)  cos(k  y  sin0  sec20)d0  [30] 


Hence,  if  we  make 


b 


n+1 


(-l)n+1kon+1(2n+l)! 
(n+1 ) ! k  2n+3 


a2n+l 


[31] 


the  wave  due  to  the  odd  power  term  of  the  source  distribution 
will  be  completely  cancelled.  The  above  argument  holds  for  any 
integer  n.  The  stern  waves  due  to  the  terms  of  odd  power  in 
the  source  distribution  can  be  dealt  with  exactly  in  the  same 
manner  as  those  of  the  bow. 
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In  general,  the  bow  waves  have  stronger  sine  components 
than  cosine  components,  the  former  depending  mostly  on  the  en¬ 
trance  angle  of  the  ship  bow.  Therefore,  the  doublet  distri¬ 
bution  is  more  important  than  that  of  the  quadrupole.  In  fact, 
this  is  the  reason,  which  will  be  shown  later,  that  the  combination 
in  practice  of  optimum  distributions  of  doublets  and  quadrupoles 
can  produce  a  closed  body,  although  the  quadrupole  Itself  can  not. 

NO-WAVE-SINGULARITY 

By  Krein's  proof  (1955)  we  know  that  there  exists  no 
singularity  distribution  which  produces  a  finite  closed  body  and 
yet  produces  no  wave.  However,  it  would  be  very  useful  if  we 
knew  of  some  singularity  which  by  itself  may  not  produce  a  closed 
body,  but  may  deform  some  closed  body  produced  by  another  singu¬ 
larity  distribution  when  both  are  combined  together,  and  which 
in  combining  produces  no  additional  wave  resistance.  Krein  (1955) 
first  found  such  a  special  function  h(x,z)whlch  consists  of  an 
arbitrary  nice  function  $(x,z)  of  x,z. 


h(x,z) 


[32] 


with  boundary  conditions 


[3^1 


on  the  boundary  of  the  domain  of  the  singularity  distribution  on 
the  center  plane  of  a  ship.  He  proved  that  the  Michell's  wave 
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resistance  due  to  a  ship  surface  f(x, z)  is  exactly  the  same  as 
that  due  to  a  ship  surface  f (x,z)  +  h(x,z).  That  is,  there  are 
infinitely  many  ship  body  forms  which  have  the  same  resistance 
and  volume. 

Here  we  will  see  that  if  we  use  a  point  quadrupole  In  the 
x  direction  and  a  point  doublet  in  the  z  direction  at  the  same 
point,  the  combined  singularity  will  not  produce  any  wave  in  far 
downstream.  The  wave  height  due  to  a  doublet  in  the  z  direction 

at  the  point  (0,0, z  )  is 

1 

tt/2 

=  8k®  J  p.  exp(-kQZ  sec20)sec50  c.os(k  x  sec  0)  X 
o 

X  cos  (k  y  sin  9  sec2©)  d0  [351 

i 

Hence  if  we  compare  Equation  [25]  with  Equation  [35]  it  will  be 
easily  seen  that 


C  .  +  C  =  0 
-d  sq 


[36] 


for 


k  X  =  p 
o 


[37] 


Since  this  is  a  point  singularity  and  it  does  not  produce  any 
wave,  it  can  be  distributed  in  an  arbitrary  manner  anywhere  on 
top  of  the  ship  singularity  distribution  without  changing  the 
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original  ship  wave  resistance. 

In  fact,  it  can  be  shown  that  if  this  no-wave-singularity 
is  distributed  continuously,  this  will  be  in  effect  the  same  as 
Krein's  form.  Besides,  by  differentiating  both  wave  heights 
and  £  by  the  same  parameter  we  may  be  able  to  construct  in¬ 
numerable  higher  order  no-wave-singularities.  Therefore  in  this 
report  we  will  call  the  aforesaid  no-wave-singularity  as  that  of 
the  first  order. 


BODY  STREAMLINE  SHAPE  DUE  TO  COMBINATION  OF  SINGULARITIES 
IN  THE  UNIFORM  FLOW 


Now  to  utilize  the  quadrupole  in  practice  we  have  to  in¬ 
vestigate  the  way  in  which  it  produces  bodies  when  combined  with 
other  singularities  in  the  uniform  stream. 

To  be  more  general,  we  may  further  have  to  consider  the 
discontinuity  (shoulder)  in  the  tangent  of  the  waterlines.  Then 
not  only  the  negative  doublet  but  also  the  positive  doublet  at 
the  slope  discontinuity  may  be  needed  to  cancel  the  negative  sine 
waves . 

In  the  following  subsections  we  will  consider  body  shapes 
which  may  be  produced  by  combinations  of  various  singularities 
which  may  be  utilized  in  improving  ship  shapes  in  order  to  de¬ 
crease  their  wave  resistance. 

The  streamlines  produced  by  the  three  dimensional  singularities 
in  the  uniform  flow  may  be  obtained  by  solving  the  streamline 
equations 


dx  dy  dz 
u  ~  v  ~  w 


[38] 
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where 


u 


Sx 


V,  v 


c)<t>  c)<t> 

ay'  w  _  "  az 


[39] 


represented  by  distributions  of  singularities.  This  could  be 
solved  by  Runge-Kutta-Gills  (Gill,  1951)  numerical  method. 
However,  the  exact  form  of  0  due  to  the  singularity  under  the 
free  surface  even  under  the  assumption  of  the  linear  boundary 
conditions  Is  so  complicated,  that  it  becomes  practically  very 
advantageous  to  relate  the  exact  shape  of  the  ship  to  the  source 
distribution,  on  the  assumption  of  zero  Proude  number.  For 
example,  Inui  (1957)  used  this  approximation  (double  model);  he 
verified  the  sufficient  validity  of  his  method  for  moderate 
Proude  numbers  by  comparing  theory  and  his  experiments.  This 
method  could  be  used  here  to  find  the  relation  between  singu¬ 
larity  distributions  and  ship  shapes,  as  in  the  case  of  ships 
with  no-wave  singularities  and  with  bulbs  made  of  doublets  and 
quadrupoles.  However,  considerable  Insight  is  gained  with  much 
greater  simplicity  by  considering  only  special  simple  cases,  as 
the  two-dimensional  case  for  instance.  Herein  we  discuss  the 
problem  using  dimensional  quantities  unless  specified. 

TWO-DIMENSIONAL  DOUBLET  AND  QUADRUPOLE 

The  complex  potential  due  to  a  doublet  with  the  strength 
-  u  <  0  and  a  quadrupole  with  its  strength  A  >  0,  with  both 
directions  parallel  to  the  x  axis,  in  a  uniform  flow  is 
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w  =  <t>  +  lY=  -  Vz  -  — 

Z  ,2 


[40] 


where 


z  =  x  +  iy 


Hence  the  stream  function  is 


Y  =  -  Vy  + 


J ±2L 


2Axy 


2  2 
x  +y 


(x2+y2  )2 


=  --Y . - 

x2  +y2 


V  (x2+  y2)  -  u  + 


2Ax 

2  2 

x4Sf+  y 


[41] 


If  we  use  polar  coordinates 

x  =  cos  0,  y  =  •  sin  0 


we  can  easily  see  that  the  dividing  streamline  is 

y  =  0 


and 


2 


jl  A.  cos  6 

V  "  V  r 


[42] 


This  means  that  the  body  is  deformed  from  the  sphere 


[43] 
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extending  (or  shrinking)  r  by  an  approximate  amount  of  —  cos  9 

in  each  direction.  In  order  to  have  a  closed  body  such  that  all 
streamlines  inside  the  surface  streamline  are  wholly  contained, 
the  magnitude  of  X  should  not  be  too  large.  To  obtain  the 
limiting  value  of  X,  we  have  only  to  consider  the  point  0=0. 

We  put 


f(r) 


[44] 


and  notice 


f(r)  =  2  ^  >  0  when  both  y  =  C  and 


And 


f '  (r )  =  0  at 


Hence  if 


f 


+ 


2 X 
V 


<  C 


or 


[45] 


there  exist  two  roots  of  f(r) 


0  in  0  < 


Tha  t  is. 


when  X  becomes  large  from  X  =  C  within  the  limit  of  the  inequality 
(45)>  Equation  [44]  represents  two  separate  closed  streamlines 
such  that  one  is  inside  the  other  (see  Figure  5a).  If  X  becomes 
larger  than  (  44-  ,  the  inner  streamlines  which  were  inside 
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the  inner  closed  streamline  when  y  <  |  37  5  pop  out  and  there 

will  be  only  one  closed  streamline  (see  Figure  5b).  Therefore 
7\  should  satisfy  the  inequality  (45)  in  order  that  a  meaningful 
body  exist.  Figures  6  and  7  show  examples  of  the  dividing 
streamlines.  When  \  is  negative  the  body  form  will  be  reversed 
along  the  x  axis,  and  the  condition  on  the  magnitude  of  A  and  p 
is  exactly  the  same  as  in  the  previous  case. 

A  SIMPLE  SOURCE  AND  A  DOUBLET  IN  A  UNIFORM  STREAM 

Consider  a  point  doublet  with  the  strength  p  >  0  at 
z  =  x  +  ly  =  0  combined  with  a  point  source  with  its  strength 

m>o  at  z  =  zQ.  Then  the  complex  potential  w  will  be  written 

w  =  -  Vz  -  m  log  ( z - ZQ )  +  ^  [46] 


The  stream  function  is 


f  =  -  Vy  -  m  tan 


-x 

x-x 


J±2L 


2  2 
x2+  y 


[47] 


By 


x  =  r  cos  0,  y  =  sin  9 


Y  =  -  Vr  sin  0  -  m0^-  ^  S*n  6  [4t>] 

where 

I  y'M 

0  =  arc  tan  - 

1  \  x‘Xo  I 
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Hence  the  body  streamline  is  obtained  by  putting  Y  =-mir.  Non- 
dimensionalizing  Equation  [46]  by  Vh  e  if 

r2  sin  9  +  r  (0  /m-l)  +  p  sin  0  =  0  [49] 

i 

The  streamline  due  to  a  source  in  the  uniform  flow  is  well  known 
(see  e.  g.  Mllne-Thomson,  1955  P*  159).  The  combination  of  source 
and  a  positive  doublet  produces  a  neck.  The  strength  of  the 
doublet  should  satisfy 


U  -  ^8  <4  [50] 

in  order  for  the  body  streamline  to  make  physical  sense.  The 
dividing  streamline  is  plotted  in  Figure  6  as  an  example,  for  a 
special  value  of  p. 

A  SIMPLE  SOURCE  AND  A  NO-WAVE-SINGULARITY  OF  THE  FIRST  ORDER 
IN  A  UNIFORM  STREAM 

We  consider  the  no-wave-singularity,  a  quadruple  in  the 
x  direction  with  strength  A  and  a  point  doublet  in  the  z  di¬ 
rection  with  the  strength  p,  both  located  at  the  origin,  and  a 
point  source  at  zq  in  two  dimensions.  The  complex  potential 
due  to  these  singularities  is 

w  =  <J>  +  i*  =  -  Vz  +  i  ^7  +  — ■  -  m  log  (s-z^)  [51] 

Hence  the  stream  function  is 
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px  _  2Axy 

2  2  /  2  p  \  2 

X  +  y  (x‘i+  y**)- 


m  tan 


-i  y-yo 

x-x 


[52] 


Since  the  no-wave-singularity  does  not  produce  any  extra  fluid, 
y  =  -  rmr  which  is  the  dividing  streamline  for  the  half  body  due 
to  a  simple  source  in  the  uniform  stream  will  still  be  the 
dividing  streamline  in  this  case. 

If  we  nondimensionalize  Y  by  hV  where  h  =  is  the  radius 

of  the  half  body  due  to  a  point  source,  the  dividing  streamline 
can  be  written  as 


0  —  — 

,  —  1  X  \1 

-1  =  -y  -  +  - - - 

7T  / — 2  — >  \2 

(x2+  y2  Y 

where 

y  =  y/'h,  x  =  x/h,  u  =  u/ (Vh2),  >  =  x/CVh3),  9i=tan  ~iy~yo 

x  -  xo. 

Although  it  is  very  hard  to  see  the  exact  shape  of  the  dividing 
streamline  without  plotting  x  and  y  of  Equation  [53]  it  may  be 
easy  to  see  just  how  the  original  half  body  streamline  is  dis¬ 
torted  by  the  no  wave  singularity.  The  last  term  in  Equation_[  53] 

is  positive  inside  the  circle  K  [  center  at  (C,  ^  );  radius,  ^  ) 

P  P 

and  x  <  0,  and  outside  the  circle  K  and  x  >  0;  and  negative  in¬ 
side  the  circle  K  and  x  >  C,  and  outside  K  and  x  <  0.  When 

x2  +  y2  »  1  the  distortion  of  the  original  dividing  streamline 
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is  small.  The  distorted  streamline  will  be  of  the  shape  shown 
in  Figure  9.  The  no-wave-singularity  has  from  Equation  [37]  the 
property  of 


X 


Vf.  _ 

gh 


For  L  <6h  and  F^  >0.3,  X/p.  >0.5.  Hence,  for  the  moderate  F^, 

the  circle  K  may  cut  the  upper  dividing  streamline  but  not 
the  lower  one  and  there  are  two  wiggles  in  the  upper  dividing 
streamline  and  only  one  wiggle  at  lower  sne  as  in  Figure  9. 


FINITE  DOUBLET  DISTRIBUTION 

SUBMERGED  AXISYMMBTRIC  BODY  WITH  FORE  AND  AFT  SYMMETRY 

In  order  to  determine  the  effect  of  limiting  the  doublet 
distribution  to  a  finite  depth  we  examine  the  difference  be¬ 
tween  the  expressions  for  ^  as  given  by  Equation  [2C,a]  with 

Br 

that  for  (,  of  Equation  [2G,bJ. 

B 


m/2 


exp(-k  f  sec20) 

Ox 


o 


2m.(f1-f ) 
sec  0 


m 


mi 


X  sin(k  x  sec  @)cos(k  y  sin  6  sec20)  d0  [5^, a] 

X  X 

This  is  the  uncancelled  wave  system  resulting  from  cutting  off 

the  doublet  distribution  at  depth  f  .  Thus  instead  of  each  term 

1 

of  the  series  in  +■  being  zero  there  will  be  a  remainder, 

B  oB 
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ir/2 


I  bn!  r — i 

'  =  exp(-k  f  sec20)  - - - — r-r —  \ 

Rn  J  0  1  k"-1s«c2<n-1)e  L 

o  0  m=0 


n  m  2mn ,  „  „,m 

k  sec  0 (f  -f) 

o  1 


m! 


X  sin(k  x  sec  0)cos(k  y  sin  0  seca0)  d0 

l  i 


[54, b] 


The  magnitude  of  £  can  be  easily  estimated  as 

Rn 


d  n  ;  r — * 

i  <8  exp(-k  f  )  — -  \ 

’Rn  o  x>  k  n-1 

0  m=0 


b  n!  A  k  m(f  -f)m 


m: 


e.g. 


-Ro 


'Ri 


<  8  exp(-k  f  ) 

0  i 


<  8  exp(-kQfi  b  <  I  +  kQ(f  -f ) )  etc. 


M54,c] 


Per  the  usual  range  of  Froude  numbers,  C  is  the  dominant 

'RO 

term  (see  Appendix).  This  decreases  exponentially  with  increas¬ 
ing  k  f  .  Although  for  increasing  n,  i  decreases  more  slowly 
o  i  Rn 

with  increasing  k  f  ,  its  magnitude  becomes  increasingly  smaller 

°  i 

than  .  Therefore,  even  if  the  wave  amplitude  is  not  com¬ 

pletely  cancelled  out  the  doublet  distribution  still  has  a  sub¬ 
stantially  favorable  effect  on  each  term. 
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SHIP  WITH  PORE  AND  AFT  SYMMETRY 

For  a  finite  draft  ship  whose  source  distribution  is 
independent  of  z,  Equations  [20],  [21]  and  [22, b]  give  for  the 
wave  height  at  large  x  due  to  the  infinitely  deep  doublet  dis¬ 
tribution  alone 


tt/2 


o 


exp(-kosec*e)  £ 


n!  b 


n 


k  sec  0 
n=0  ° 


sin  (k  x  sec  9) 

l 


X  cos  (k  y  sin  0  sec20)  d0  [56, a] 

i 

Equation  [56, a],  of  course,  also  represents  the  negative  of  the 
wave  system  due  to  the  ship  bow  where  b^  is  obtained  from  Equa¬ 
tion  [  21]  . 

If  we  omit  the  doublet  distribution  from  the  point 

z  =  f  >  1  to  z  =  »,  we  obtain  the  uncancelled  wave  height 

XX  X 

due  to  this  by  integrating  the  second  term  in  Equation  [2G,a] 
with  respect  to  f  from  0  to  1. 


tt/2 

00  r 

r —  /  r >  r — >  I 

=8/  exp(-k  f  sec2©))  b  n! \  — - - - - 7 - t-t 

HZ>  J  01  L "  L  (r+1 )  !  sec2  <n-r-1) 


n  (f  r+1-(r  -i)r+1)k  r*n+1 


n=G  0 


n=C  r=0 


X  sin  (k  x  sec  0)cos(k  y  sin  9  sec20)  d0  [56, b] 

x  x 
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This  is  formally  similar  to  Equation  [24, a].  As  in  Equa¬ 
tions  [24,  b]  and  [24,  c],  we  can  write 


\ 


l-Ro 

<  6 

exp ( -k  f  )  b  k 

v  0  1  *  00 

L'Ri 

00 

V 

exp ( -i:  f  )  b 

O  1  1  ] 

1  i-2f. ) 

i1-  2  \ 

k 

0 

etc . 

\  1 

i 

Therefore  the  same  comment  about  Equation  [54, a]  can  be  applied 

to  Equation  [56, a]  that  is,  the  uncancelled  wave  height,  f  of 

'R 

Equation  [56,b]  decreases  almost  exponentially' with  increasing 
kQf  for  the  usual  range  of  Froude  numbers. 

To  obtain  the  approximate  correction  we  may  use  a  per¬ 
turbation  method.  We  replace  b^  by  b  (l+€  )  in  Equations  [56, a] 
and  [56, b].  The  sum  of  these  two  newly  formed  equations  gives 
the  wave  height  produced  by  the  perturbed  doublet  distribution 
of  depth  f.  When  we  subtract  from  this  the  wave  system  due  to 
the  ship  bow  as  given  by  the  negative  of  Equation  [56, a]  we  ob¬ 
tain  for  each  term  of  the  resulting  series 


I 

n 


tt/2 


o 


exp(-l<osec2e  ) ) 


seczn  0 


exp(-k  f  sec20) 
e  — ( 1+ e  ) - — - 


,ec2n  9 


^  [f  r+1-(f  -l)r+1 
V  Li _ 1 _ J 


‘E 

r=0 


r+1 


(r+1 )  I  cos2 ^0 


sin(k  x  sec0)cos(k  ysin0sec20  )d0 . 

1  1 

[57, a] 
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The  wave  resistance  represented  by  the  first  term  in  the  wave 
system  I  (putting  n  =  0  in  Equation  [57, a]  except  for  a  con¬ 
stant  factor  is 


m/2 

r  ; 

1  1 

^  -kosec20^ 

-k0f  sec2© 

1 

e 

1  ° 

l1'*  J 

- (l+e  )e  k  sec2©  > 

'o'  o 

✓ 

o 

1 

i  o  -k_  -k0/2  k 

-  e2  <  -  +  2  e  G  (k  )-4e  G  -§■  I  \  -2e  (l+e  )  k 

o)3  -a  o'  -2  \  2  \  (  o'  o'  o 


-k  f  /a  /k  f  -(k  +k  f  )/s  I  k  +k  f 

X  <e  “V,  U*  -e  ° 


+  k  2  (l+e  ) 
o'  o' 


a  e 


-k  f 
o  l 


- —  K  (k  f  ) 
2  o  o  i 


where 


t/2 


-k. 


exp(-2kQsec20  )sec2n+10d© 


G  (k  ) (See  Yim  1962a) 


,n+l  n  o 


-l 


k  f 
o  1 

\  2  / 


k  f 

o  1 


k  f 

o  1 


\  2  / 


-  K 


k  f 
o  1 


o  2 


k 


G  ,  (k  )  =  K  (k  )  +  2  k  K  (k  )  -  2k  K  (k  ) 

-2  'o'  3  \  1  o'  o  o'  o'  o  lO'< 


K  (k  ),  K  (k  )  are  the  modified  Bessel  functions  of  the 
o'  o'  i  '  o' 


second  kind  of  zeroth  and  first  order  respectively.  Hence  the 
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value  of  c  which  makes  R(I  )  minimum  is  e  =  ^ 
o  '  o'  o  2a 


where 


_  -k  -k  /2 

ot  =  —  +  cto  tt  (k  )  —  *r  e  G 

3  -2  0  -2 


-2k 


j3  =  -  2k  ( 
^  o  ' 


J  -VA  | 

1  Vx  ' 

-(k  +k  f  )/2 

°  0  2."  r. 

f  k  +k  f  \ 

Y  0.x  ! 

1 2  1 

-e  G 

l  2  j 

v  +k2  e  -  K  (k  f  ) 

•  o  2  o  o  i 


-k  f  / 2 
e  °  1  G 


k  f 
o  1 


-(k  +k  f  )/2  k  +k  f 

O  O  i  o  O  1 

-e  1  G  - - — ~ 

~1  ^  j 


-k  f 

+  k  2  e  °  1  K  (k  f  ). 
O  0  o  1 


The  values  of  e  for  different  Froude  numbers  and  f  =1  are 
o  i 

given  in  Table  1.  The  values  of  e  for  1  can  be  calculated 

n 

as  above.  However,  since  these  are  not  so  Important  in  the 
usual  range  of  Froude  number,  a  rough  and  easily  obtained  ap¬ 
proximation  may  be  sufficient. 

Using  the  method  of  stationary  phase  for  y  =  0 


I 


n 


=-8 


-k  f 
o 


e  -(1+e  ) 
n  v  n '  k 


0  L—i 


r+1-(f  -i)r+1  k  r+1 
1  ’  0  . 


(r+1)! 


X  sin  (k  x  +  77  )• 
i  + 


[57, b] 
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If  we  put  I  =  0 
^  n 


e 

n 


[57*  c] 


Equation 

numbers. 


[57>c]  is  evaluated  for  f  =1  and  different  Froude 

i 

from  n  =  1  to  n  =  3  are  shown  in  Table  1.  For  f  = 

i 


1. 


n-f-1  ,  u 

k 

o 
m! 


m=l 


m=l 


[57, d] 


We  note  that  e  in  [57, c]  is  not  a  function  of  a  or  b 
n  n  n 

but  only  of  f  ,  k  ,  and  n,  and  is  always  positive  for  f  >  1 
10  1 

since  the  numerator  in  Equation  [ 57, c ]  is  positive,  and  the  de¬ 
nominator  in  Equation  [57,c]  is  larger  than  that  of  Equation 
[57>d].  The  latter  is  positive  due  to  the  fact  that 
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1  + 


n+1 


y 

L. 


k 

is  a  partial  sum  of  the  series  expansion  of  e  °.  (This  is  also 

true  for  y  /  0  in  Equation  [ 57, a ] .  Hence  if  the  b  1 s  are  all 

n 

positive  (as  for  the  sine  ship  treated  in  the  next  Section),  each 
horizontal  section  of  the  bulb  represented  by  Equations  [21]  and 
[23,b]  in  0  <  z^  <  1  will  be  smaller  than  any  corrected  bulb. 

Since  we  determined  from  the  stationary  phase  at  y  =  0, 

en's  (n  £  1)  in  Table  1  are  overestimated,  or  rather  Equations 

[21]  and  [23*b]  with  b  (1+e  )  instead  of  b  would  indicate  the 

n  n  n 

upper  bound  of  the  size  of  bulb  for  minimum  wave  resistance  sub¬ 
ject  to  the  conditions  of  this  method. 

In  order  to  find  the  exact  optimum  doublet  distribution  on 

0  <  z  <  f  for  the  minimization  of  the  wave  resistance  due  to  a 
i 

given  ship  it  is  better  to  attack  this  problem  by  a  somewhat  dif¬ 
ferent  approach.  That  is,  we  determine  the  unknown  coefficients 
b^  of  the  doublet  distribution  directly  by  minimizing  the  wave 
resistance  of  a  given  ship  with  the  bulb  (see  Yim,  1962b).  This 
is  equivalent  to  the  present  method  provided  we  obtain  en  by 

minimizing  the  wave  resistance  due  to  the  complete  uncancelled 
wave  system  directly.  is  not  necessarily  always  positive, 

since  in  this  case  there  are  many  extra  cross  product  terms. 


38 


HYDRONAUTICS,  Incorporated 


According  to  Figures  la  and  2a,  the  exact  optimum  bulb  is  thinner 
near  the  free  surface  and  thicker  near  the  keel  than  the  one  cor¬ 
rected  by  the  method  given  in  this  paper. 

It  should  be  mentioned  here  that  there  is  no  interaction 
between  the  bow  and  the  stern  waves  for  the  ideal  doublet  dis¬ 
tribution  since  in  this  case  the  bow  and  stern  waves  are  separately 
cancelled  out  by  the  bow  and  stern  bulbs  respectively.  However 
for  the  finite  doublet  lines  there  exist  interactions  and  the 
curves  of  wave  resistance  versus  Froude  number  have  humps  and 
hollows.  Therefore  the  design  Froude  number  of  the  ship  would 
normally  be  selected  so  that  it  falls  at  a  hollow  on  the  resist¬ 
ance  curve.  However,  since  the  optimum  bulb  at  the  bow  or  the 
stern  for  each  Froude  number  has  the  effect  of  smoothing  out  the 
resistance  curve  to  a  great  extent  (Ylm,  1962a)  the  unfavorable 
effect  of  not  falling  on  a  hollow  is  usually  outweighted  by  the 
advantage  of  using  an  optimum  bulb. 

CASE  OF  SINE  SHIP 

Inui  (i960)  pointed  out  that  the  smaller  the  Froude  number, 
the  greater  becomes  the  Importance  of  the  first  term  in  Equa¬ 
tion  [16].  In  the  practical  case  we  need  not  take  many  terms  in 
the  series.  At  moderate  ship  speeds  even  one  term  will  be  enough 
since  the  second  term  of  the  series  [lb]  produces  wave  heights  of 

the  order  1/k  2  <  0.01  for  F  =  .3  while  the  first  term  is  con- 
1 

sidered  to  be  of  order  1  (see  Appendix).  This  is  natural  since 
it  is  well  known  that  the  first  term  of  the  series  [16]  for  the 
source  distribution  of  the  ship  is  proportional  to  the  angle  of 
entrance  of  a  ship  and  the  angle  of  entrance  has  a  great  influence 
on  the  wave  resistance. 
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As  an  example,  the  first  four  terms  of  the  cosine  series 
will  be  taken  for  the  source  distribution  in  Equation  [16], 

Then  the  corresponding  ship  will  be  a  sine  ship  in  the  Michell's 
sense.  The  coefficients  a^  are  taken  such  that  the  ship  half¬ 
breadths  are  nondimenslonalized  with  respect  to  the  half  beam  and 
the  stations  are  measured  from  the  bow  and  nondimensionalized  with 
respect  to  the  ship  length  L.  The  values  of  a^,  b^  from  Equa¬ 
tion  [21],  and  from  Equation  [57, c]  for  different  Proude  num¬ 
bers  and  length-draft  ratios  are  shown  in  Table  1.  The  corres¬ 
ponding  doublet  distributions.  Equations  [21]  and  [23>b]  without 
perturbation  corrections,  in  0  <  :  <  1,  for  different  Proude 

i 

numbers  are  plotted  as  the  solid  curves  in  Figures  la  and  2a. 

Also  shown  in  these  figures  for  purpose  of  comparison,  are  cor¬ 
responding  curves  determined  from  the  exact  ship  wave  resistance 
minimization  theory  based  on  finding  directly  the  optimum  doublet 
distribution  for  a  given  ship  with  a  bulb  which  does  not  extend 
below  the  keel  (see  Yim,  19o2b). 

The  two  dimensional  approximate  relation  of  bulb  radius  r 
and  the  strength  of  doublet  u  per  unit  length  of  z  is  q  =  r2V. 

This  relation,  although  overestimating  the  actual  radius  r,  is 
used  as  an  alternate  means  of  expressing  the  nondimensional  doublet 
strength  in  Figures  la  and  2a  in  order  to  give  some  idea  of  the 
approximate  area  distribution  of  the  bulb. 

The  remarkable  reductions  of  wave  resistance  due  to  the 
bulbs  for  both  cases  are  shown  in  Figures  lb  and  2b.  The  bulb 
sizes  near  the  keel  (z  =  -l)  determined  from  Equation  [21]  are 
very  much  smaller  than  those  shown  by  the  exact  method  illustrated 
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by  the  broken  curves  although  near  the  free  surface  the  former 

are  slightly  larger  than  the  latter.  The  Proude  number  of  least 

resistance,  where  the  reduction  of  wave  resistance  is  almost  90 

percent,  is  usually  less  than  the  Froude  number  for  which  the 

calculation  of  the  doublet  distribution  [21]  was  made.  This  is 

due  to  the  fact  that  we  have  cut  off  the  infinite  ideal  doublet 

distribution  at  the  keel.  The  figures  of  doublet  distributions 

[21],  uncorrected  by  perturbations,  appear  to  be  almost  linear 

in  0  >  :  >  -  1.  This  indicates  the  importance  of  the  first  term 

b  a  of  the  series  for  the  doublet  distribution  [21]  or  the  first 
o  l 

term  a  of  the  series  for  the  source  distribution  for  a  ship  in 
o 

the  practical  range  of  Froude  numbers,  as  pointed  out  earlier  and 
in  the  Appendix.  That  is,  the  equation  for  the  first  term  in  the 
doublet  series  bQ  =  aQ/k  is  a  very  simple  but  very  important 

approximate  relation  between  the  skip  and  the  corresponding  bulb. 
It  is  clear  from  the  equation  for  the  perturbation  [5?]»  and 
Table  1  that  the  nonlinearity  Increases  when  the  length  of  bulb 
becomes  smaller.  Hence  it  may  be  desirable  to  increase  the  mag¬ 
nitude  of  the  doublet  strength  near  the  keel  or  rather  extend 
the  doublet  distribution  below  the  keel  if  possible  according  to 
the  methods  indicated  in  Equation  [ 5T 1 •  However,  according  to 
Figures  lb  and  2b  the  doublet  distribution  of  Equation  [21]  which 
gives  a  much  smaller  bulb  than  the  exact  optimum  one  results  in 
a  very  favorable  effect  even  without  the  perturbation  corrections. 
Besides  even  these  bulbs  are  quite  large  when  the  Froude  number 
is  large.  In  fact  we  may  expect  that  the  effects  of  viscosity. 
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cavitation,  and  other  design  conditions  would  ultimately  deter¬ 
mine  the  maximum  bulb  size  at  these  Proude  numbers. 

CASE  OP  THE  PARABOLIC  SHIP 
When  we  consider  the  source  distribution 


m  =  c  (1  -  2x)  in  0  <  z  <  1 

i 

Michell's  linear  ship  relation 


m(x) 


1  5y 
27 r  dx 


will  indicate  that  this  is  a  parabolic  ship  with 


[58] 


c 


1  B 
v  L 


where  B/L  is  the  beam-? ength  ratio. 

Since  we  have  an  even  power  and  an  od  •  power  term  in  Equa¬ 
tion  [58],  the  bow  waves  contain  both  the  '.ne  and  cosine  com¬ 
ponents.  To  cancel  the  sine  bow  waves  we  need  to  put  a  doublet 
line  of  the  strength  from  Equation  [23, b] 


£ 

c 


0  <  z  <  1 
i 


i 


1  <  z 

1 
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The  effect  of  doublet  lines,  both  the  ideal  distribution  cut  off 
at  the  keel  and  the  optimum  finite  distribution  has  already  been 
demonstrated  for  a  sine  ship.  In  Figures  3a  and  3b  the  optimum 
finite  doublet  distribution  for  the  parabolic  ship  and  the  cor¬ 
responding  bow  wave  resistance  with  and  without  the  doublet  dis¬ 
tribution  are  shown  (the  process  of  calculation  is  discussed  in 
Yims  paper  of  19o2b). 

To  cancel  the  cosine  bow  waves  we  need  to  superimpose  the 
quadrupole  line  of  the  strength  from  Equations  [28]  and  [29], 


X 

c 


in  0  <  z  <  1 
i 


X 

c 


(2  -D21 
1 


(2z  -1) 

l 


i 


in  1  S  z 

i 


k  is  usually  larger  than  1  and  quite  large  if  Froude  number  is 
a. 

small.  is  smaller  by  the  factor  of  the  draft  length  ratio 

than  k  .  Hence  inequality  (45)  would  be  satisfied  for  quite  a 
i 

large  z  for  moderate  Froude  numbers  but  eventually  inequality 
i 

(45)  is  violated  for  z  larger  than  a  certain  number.  However 

i 

the  effect  of  cutting  the  line  off  at  a  large  z  may  not  be 

i 

serious  as  was  shown  for  the  case  of  doublet. 

The  wave  resistance  due  to  the  cosine  bow  waves  of  the  para¬ 
bolic  ship  with  and  without  the  ideal  quadrupole  line  cut  off  at 
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the  keel  are  shown  In  Figure  4b.  The  optimum  quadrupole  strength 
and  the  wave  resistance  with  and  without  the  optimum  quadrupole 
(see  Yim  1962b)  is  also  shown  in  Figures  4a  and  4b.  Although  the 
effect  of  the  quadrupole  line  is  not  as  prominent  as  that  of  the 
doublet  for  the  sine  wave,  even  the  ideal  quadrupole  line  cut  off 
at  keel  is  shown  to  be  quite  favorable  for  moderate  Froude  numbers. 
The  wave  resistance  due  to  the  cosine  component  of  the  bow  waves 
even  without  the  bulb  as  is  shown  in  Figure  4b,  is  quite  small 
compared  to  that  due  to  the  sine  component.  However,  in  practice 
the  cosine  components  seem  to  be  much  larger  than  in  this  theory 
due  to  the  change  of  the  actual  body  streamlines  caused  by  boundary 
layers  and  wakes. 
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APPENDIX 


Since  a  is  the  coefficient  of  a  power  series  which  re- 
20 

presents  the  source  distribution  of  a  ship,  it  is  usually  com- 

2n 


parable  to  the  cosine  series  where  a 

2n 


(-if 


7 r _ 

(2n)! 


It 


therefore  seems  reasonable  to  assume  that  (2n)l  a 

2n 


<  Mtt2 


n 


M  = 


oil 

finite  number,  for  all  n.  For  usual  ships  F  =  —  <  .  Be¬ 
sides  the  Integral  1 


C  = 


tt/2 

/ 


cos 


2n-i 


0d0 


2.4.  - .  2  (n-1 ) 

3.5.  —  .  (2n-l) 


is  decreasing  with  increasing  n.  Hence  a  typical  term  in  the 
Equation  [l8]  for  the  wave  height  due  to  the  ship  can  be  esti¬ 
mated  as 


tt/2  -k  f  sec20  (2n)la 

/  e  °  - 2-9 — — »  sin(k  x  sec  0)cos(k  y  sln0sec20)d0 

o  (k  sec  ‘  1 

1 


_k  f  _  V2  „^„2  (n-1 ) 


2  -k  f 

.x  „  2  ,  cos-'  '0  r  MC  .  o  r.,, 

O  Mtt  /  .  /„  i  \  d0  -  — 7T  i  \  e  [Al] 


<e  "  o'  (k/ir)2(n-l)  '  ^2  (n-1) 


(k  /rrY 


Since  7r/k  is  a  small  number  usually,  it  is  clear  how  rapidly 
i 

the  integral  [Al]  decreases  when  n  increases.  Hence  in  practice 
we  do  not  need  to  take  many  terms  of  the  series  in  determining 
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the  optimum  bulb.  That  is,  the  magnitude  of  the  corresponding 
terms  for  the  doublet  distribution  decrease  very  rapidly  in  the 
usual  range  of  Proude  number.  Therefore  we  have  only  to  con¬ 
sider  the  case  of  finite  n  when  we  consider  the  limiting  case  of 
2  _►  co  •  i.e.  the  effect  of  the  infinitely  long  doublet  distri- 

i 

button  from  Equation  [20, a].  For  the  case  of  a  finite  doublet 
distribution  we  have 


-k  f 
o 


.  n-1 
k  sec 
o 


sec20 
2 (n-1  )9 


-k  f1seca6 
o 

-e 


I 


(f  -f)n-r 

(n-r)ik  r"1sec2^r"1^0 
o 


sin(k  x  sec  0) 

i 


X  cos  (k  y  sin  0  sec20)  d9 
i 

Let  us  denote  the  terms  due  to  the  limit  f  by  R  . 

i  n 

finite  n 

lim 


f  —  00  R  —  0 

a  n 


Then  for  a 


because 


11m  -k  f 

o 

f  -  c»  e  (f1-f)n=o 


The  argument  is  the  same  when  the  amplitude  functions  are  given 
by  Equation  [ 22 ] . 
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NON-DIMENSIONAL  D'T'JBlET  STRENGTH,  i2»AhBV), OR  SECTIONAL  AREA  I2r2/(BM> 

FIGURE  20-DOUBLET  DISTRIBUTION  FOR  SINE  SHIP  (L/H  =  24) 


FROUDE  NUMBER  FL 


FIGURE  2b-  BOW  WAVE  RESISTANCE  OF  SINE  SHIP  (  L/H  =  24) 
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FROUDE  NUMBER 

FIGURE  3b  BOW  WAVE  RESISTANCE  DUE  TO  THE  EVEN  TERM  IN  THE 
SOURCE  DISTRIBUTION  OF  A  PARABOLIC  SHIP 
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NONDIMENSION  AL  QUADRUPOLE  STRENGTH, 

(VBH  ) 


FIGURE  4a  QUADRUPOLE  DISTRIBUTION  FOR  A  PARABOLIC  SHIP,  (-^  =  16) 
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FROUOE  NUMBER 


»b  BOW  WAVE  RESISTANCE  DUE  TO  THE  ODD  TERM  IN  THE 
SOURCE  DISTRIBUTION  OF  A  PARABOLIC  SHIP,  ( =16) 
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FIGURE  5  SCHEMATIC  DIAGRAM  OF  DIVIDING  STREAMLINES  DUE  TO 
A  POINT  DOUBLET  PLUS  A  POINT  QUADRUPOLE  AT  THE 
ORIGIN  WITH  THE  DIRECTIONS  PARALLEL  TO  THE  AXIS 
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FIGURE  8  TWO  DIMENSIONAL  DIVIDING  STREAMLINE  DUE 
TO  A  POINT  SOURCE  AND  A  POINT  DOUBLET 
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